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Abstract 

We calculate the molecular-spin (5), temperature (T), and field (H) dependence of 1/T% for a 
local magnetic probe coupled to an isotropic high spin molecule, based on spin-phonon interaction. 
We compare the calculation to recent NMR and /iSR experiments in CrCu6 (S = 9/2), CrNig 
(S = 15/2) and CrMng (S = 27/2). Although we can account for the high and intermediate 
temperature regimes, the calculation is fundamentally different from the data at T — ► 0. Since 
1/Ti must be due to coupling of the molecular spin to an external heat bath, and since phonon 
contribution is ruled out at low T, we conclude that at these temperatures hyperfine interactions 
must play an important role in the molecular spin dynamic. 
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I. INTRODUCTION 



High spin molecules (HSM) consist of clusters of metal ions, they are ordered in a crystal 
lattice, and coupled by Heisenberg ferromagnetic or anti-ferromagnetic interactions with 
coupling constant J, only between spins Si in the molecule. In this paper we study a family 
of isotropic high spin molecules. We calculate the spin lattice relaxation 1/Xi in these 
molecules, and compare it to experimental measurements^ to determine the origin of the 
observed spin dynamics of the molecules. We find that at high temperature the molecular 
spin dynamics is driven by thermal activation (spin-phonon interaction). However at low 
temperatures the molecular spin dynamics is dominated by hyperfine interactions between 
the molecular and nuclear spin. 

These calculations can be easily applied to different kinds of high spin molecules, such as 
Mn 12 i'iS and Fe§i, where quantum tunneling of the magnetization (QTM) is observed&Si. 
Comparison between experimental measurement si'B and the calculation can help determine 
which interactions induce the spin tunneling in these molecules. 

In the isotropic molecules a Cr(III) ion is surrounded by six cyanide ions, each bonded to a 
Cu(II), Ni(II) or Mn(II) ion. The coordination sphere of Cr and Cu/Ni/Mn can be described 
as a slightly distorted octahedral. For convenience we will refer to these molecules as CrCu6, 
CrNi 6 and CrMn 6 , respectively. The Hamiltonian of the isotropic high spin molecules at zero 
field can be written as 

h = -jEvs« (1) 

i=l 

where Sq is the spin of the central Cr ion (with S = 3/2), i runs over the peripheral Cu, Ni 
or Mn ions (with 5* = 1/2, S = 1 or S = 5/2), which are coupled to the central Cr ion, with 
coupling constant J (e.g. see Figure |I|). At temperatures lower than J only the ground spin 
state S = 9/2, S = 15/2 or S = 27/2 is populated. 

The Hamiltonian 7io is isotropic therefore the total spin S and the spin in the z direction 
m are good quantum numbers, and the eigenfunctions of Tio can be written as \S,m >. 
However, at very low temperatures T < J, the degeneracy of the ground state is removed 
by additional anisotropic perturbation on 7io- Such perturbation that does not commute 
with Tio can cause transitions between the different spin states \S, m >aa 
and induce the observed spin dynamicsEE. 

The anisotropic term in the Hamiltonian can be written as the sum 7i l = H c + H n , where 




FIG. 1: The magnetic core of the CrNig molecule, So is the spin of the Cr ion and Si, S2, • • • , Sq 
are the spins of the Ni ions. 

7i c commutes with the Hamiltonian Hq, while 7i n does not. These terms may be a result 
of dipolar interaction between neighboring spinsHH, spin-phonon interactionEl'0, nuclear 
fluctuationsEl, high order crystal field termslli'i00'0 or small anisotropy in the coupling J 
between spinsEl. 

In order to calculate the value of the spin lattice relaxation time T\ in the isotropic HSM 
we diagonalize the Hamiltonian Tt + 7i c in Section [TJ. In Section [TIJ we show that we can 
calculate the magnetization and susceptibility of the compounds using the eigenvalues and 
eigenfunctions that we obtain. In Section we use the eigenvalues and eigenf unctions of 
the Hamiltonian to calculate 7\, taking 7i n as a perturbation. The perturbation introduces 
a finite lifetime for the different spin states \S, m > and induces spin dynamics, or transitions 
between the different spin states. 

II. EXACT DIAGONALIZATION OF THE HAMILTONIAN 

To calculate the eigenvalues of the Hamiltonian (HD we write it in the form 
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H = -JS -^2 ( 2 ) 

i=l 

Using S = S + Yh=i Si and St = J2t=i $i one can WT ^ e 

Ho = -JSo -S t = - J - (5 2 - 5? - Si) . (3) 
The energy eigenvalues of the sates \S,m >= \So, S t , S,m > are 

E s = E(S, S t , So) = - 3 - (S(S + 1) - S t (S t + 1) - So(S + 1)) (4) 



where we use the notation S for the set of numbers (S, S t , So). The degeneracy of the state 
\S, m > is the degeneracy of the value of St- When an external magnetic field H is applied 
the Zeeman term should be added to 7Yo, and the full Hamiltonian becomes 



Tt = Hq — gfisHSz, 
for which the eigenfunctions \S,m > are not changed and the eigenvalues are 
E s , m = ~ (S(S + 1) - S t (S t + 1) - S (So + 1)) - gfJL B Hm. 



(5) 



(6) 



In Figure |2| we present the energy levels of CrNig at zero magnetic field H = as a 
function of the total spin S (a) and as a function of the spin in z direction m (b). 
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FIG. 2: The energy levels of the Hamiltonian Ho at zero field for CrNi6 as a function of (a) the 
total spin S and (b) the spin in the z direction m. 



III. SUSCEPTIBILITY CALCULATION 

The susceptibility can be calculated as a function of temperature and external field using 
the eigenstates and eigenvalues that we calculated in the previous section. First we calculate 
the magnetization 



TUP T 

M(T) =< S z >= £ (7) 

\S,m> 

where Z = J2\s,m> ex P ( — ""r 12 ") i s the partition function. The measured susceptibility! is 
defined as 

X(T) = (8) 



In Figure |^ we fit the experimental measurement of xT as a function of temperature in (a) 
CrCu6, (b) CrNi6 and (c) CrMng, at fields H = 100 G and 2.15 T, to the calculated value of 
\T from the Hamiltonian ([|) and Eq. (S). From the fit one obtains the coupling constants 
^o-Cu = 77 K, Jo-Ni = 24 K and J C r-Mn = — 11 K and no anisotropy is observed in all 
three molecules, within the fitting accuracy. This indicates that the high spin molecules can 
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FIG. 3: The susceptibility multiplied by temperature as a function of temperature at two different 
external fields, measured in (a) CrCug (b) CrNig and (c) CrMng. The solid lines are fits to the 
theoretical expectation (see text). 

be well described by the Hamiltonian ([5]), and that the molecules are indeed isotropic. 



IV. Ti CALCULATION 

The spin lattice relaxation time T\ was measured^ in CrCu6, CrNig and CrMng using /zSR 
(where the local probes are polarized muons) and proton-NMR. Assuming for simplicity an 
isotropic interaction between the probe I (muon or proton) and the local electronic spins 5* 
of the whole molecule 

His = AI-S (9) 

the spin lattice relaxation is given bylll 

i /l 2 r°° 

— = — (S-{t)S+(0)) e wt dt (10) 



Ti 2 J- 

where uo = is the probe's Larmor frequency in an external field H. The value of 

(S-(t)S+(0)) = (e- im/h S^e mt/n S + ) 
for the eigenstates \S, m > of the Hamiltonian TC in Eq. (^) is, 

,( g S,m+l- g S,m \ , p 

(S-(t)S+(0)) = Yl ( S ( S + 1) - m(m + 1)) e V fi J - 



\S,m> 
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FIG. 4: The values of the spin lattice relaxation measured by ^uSR (l/T-^), compared to those 
measured by proton NMR (l/Tf 1 ) after scalling, at the same external field H = 2.15. 

where Z is the partition function. This result is correct for a system with levels of infinitely 
long lifetime. However, assuming a Lorentzian broadening of the levels, due to the non- 
commuting additional terms in the Hamiltonian, 7i n , we arrive at 

E S,m 

(S-.(t)S+(0)) = ]T 6 J (S(S + l)-m(m + l))e~^e' 



\S,m> 



Z 



where Ts, m is the lifetime (or the inverse of the broadening) of the level \S,m >. Hence 
according to Eq. (|T0l) 



Y = Y E ^-{S{S+l)-m{m+l)) 

\S,m> 



(12) 



where w' = uj + (J5s, m +i — Es >m )/h. The energy difference in our case of isotropic molecules 
is -E's.m+i — Es m = —gfi B H, therefore uj' = uj — S ^-H and 



Y = ^ E (5(5 + l)-m(m + l))(I fl - me " 

' |5,»n> > 



+ W T S,m 



(13) 



in contrast to what was obtained inS, where the term (S(S + 1) — m(m + 1)) is missing in 
Eq. ( |i~3"D and uj' is replaced by uj = jH. 

This result is important for understanding of the comparison between the spin lattice 
relaxation rate measured by /iSR and NMR at H = 2.15 Tesla, shown in Figure f|. The 
field dependence of 1/Ti comes from uj' = (7 — g/j,B/fi)H. Since the gyromagnetic ratio 7 of 
the probe (muon or nucleus) is much smaller than g/xe/fo (electronic gyromagnetic ratio), 



one can write u' ~ —gfisH/h which does not depend on the value of 7. Therefore the 
field dependence of the spin lattice relaxation rate l/Ti is independent of the value of 7, 
and depends only on the value of A. This explains the fact that the spin lattice relaxation 
rate measured by proton-NMR can be scaled to match that measured by pSR, at the same 
external field, and not at the same Larmor frequency of the probe (see Figure However, 
Eq. flTB| ) is valid assuming that the perturbation 7i n is smaller than the Zeeman splitting, 

<m+ \ — Es, m , i>e. at high fields where gpsH ^> Ti n . When the Zeeman splitting is smaller 
than TC n , uj' should be replaced by u = jH, in Eq. (Dp. 

The lifetime Ts )in of the level \S, m > can be expressed in terms of transition probability 
from the state \S,m> to another state \S',m' > 

J_ = ]T p(S,m->S',m') (14) 

which depends only on the additional parts of the Hamiltonian 7i n which induce these 
transitions. In what follows we will try to account for this lifetime assuming different possible 
interactions. 



A. Spin-Phonon Interaction 

To account for the temperature dependence of T\ we should take into considera- 
tion the transitions induced by spin-phonon interactions. The spin-phonon coupling 
Hamiltoniai]0@, Tin, can induce transitions between different spin states of the molecule. 
The transition rate from a state \S, m > to a state \S',m' > can be calculated using the 
golden rule in perturbation theory^ 

fa of a 3 I (S,m\H sp \S',m') | 2 1 
p(S,m ^ S ,m) = (E Sim -E s - 1 



2vr h A pc 5 v °' m S '' m ' J ex V [(Es. m -E S i, m i)/T}-l 

(15) 

This result involves the matrix element (S, m\H sp \S' ', m') of the spin-phonon interaction, the 
phonon velocity c, the specific mass p and the energy difference (Eg >m — Es' jTn '), where it 
was assumed that Es, m > •£'s',m'- To get the right order of magnitude, and simplify the 
calculations, we assume a constant spin-phonon interaction matrix element, arriving at 

p{S 1 m^S,m) = -— w ^ = -j-^ — - (16) 
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where 

3 | (S, m\T-C sp \S', m') 



C 



2tt h A p& 

The transition probability due to spin-phonon interaction strongly depends on tempera- 
ture. At very low temperature phonons die out exponentially with decreasing temperature, 
yielding a very low transition probability, and extremely low spin lattice relaxation rate 
values as seen in Figure [5]. Therefore this interaction cannot give a full explanation to the 
nonzero spin lattice relaxation rate at low temperatures, which is observed in experiments&lll, 
and additional terms in the Hamiltonian should be considered. 




FIG. 5: Spin lattice relaxation rate as a function of temperature for different magnetic fields in 
CrNig, when assuming spin-phonon interactions only. For this figure we used C = 400 1/sec K 3 
and A = 5.2 MHz in Eq. (pi) and (0). 



B. Other Interactions 

In order to obtain a finite spin lattice relaxation rate at very low temperatures, the lifetime 
of the levels should be finite. This cannot be accounted for by spin-phonon interaction as 
seen in the previous section. A finite lifetime can be achieved if one simply assumes a finite 
broadening of the levels due to an additional interaction 7ii n t, giving a short lifetime Tj nt for 
the levels. The assumption implied by the experimental results^ is that Ti nt is temperature 
and field independent at low fields. 

In this case the total lifetime of the levels consists of two contributions, r sp due to spin- 
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phonon interaction H sp and T int due to the additional interaction Hint, therefore 



1~S,m \ Tsp Tint 

At high temperatures the value of the spin-phonon contribution to the lifetime, r sp , is much 
shorter than r int and the value of 7\ is dominated by spin-phonon induced transitions, while 
at very low temperatures r sp is much longer than r int , and the value of Ti is dominated 
by Ti n t- At intermediate temperatures both contributions to the spin lattice relaxation are 
important. 

In Figure || we present the experimental values of 1/Ti, measured by /iSR@, as a function 
of temperature at different fields, for the molecules (a) CrCu6, (b) CrNi 6 , and (c) CrMn 6 . 
The solid lines are fits to the theoretical values expected assuming transitions which are 
induced by spin-phonon interaction in addition to rii nt . The fits give the parameters' values 
summarized in Table |. 



Compound 


r int [nsec] 


C [1/sec K 3 ] 


A [MHz] 


CrCu 6 


7.0(8) 


0.13(9) 


1.7(1) 


CrNi 6 


11.0(8) 


400(60) 


5.2(2) 


CrMn 6 


9.1(8) 


0.004(1) 


4.7(2) 




TABLE I: The fit parameters of the theoretical calculation of the spin lattice relaxation to the 
experimental values from ^SR measurements 

The parameters quoted in Table | were calculated using Eq. ( p~3| ) with u = j^H instead 
of u/, since the experimental measurements were performed at low fields. The quoted values 
of T^t indicate that rii nt is of order of 0.7— 1.1 K, which is larger than the Zeeman splitting 
in fields up to 2 kG, and self consistent with the use of Eq. (|IB|) with uj = 'y fl H instead of to'. 

The fits in Figure || capture the essence of the temperature and field dependence of 
1/Ti, considering the simplifications that we have used, and it gives the correct general 
behavior of the experimental data. However, at very low temperatures and high fields the 
theoretical calculation deviates from the experimental data. Similarly at high temperatures 
T ^> J, where 1/Ti is very small, the theoretical calculation deviates from the experimental 
data (especially in the case of CrNig). We believe that the origins of the deviation at high 
temperatures is that the Hamiltonian (§) does not describe the system well enough, and 
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FIG. 6: The value of 1/T\ as a function of T at different fields, measured using /xSR in (a) CrCug, 
(b) CrNig and (c) CrMng. The solid lines are fits to the calculated value (see text). 

that the value of 1 jT\ is very small and is harder to estimate experimentally. This can also 
be seen in Figure [|, where at high temperatures the calculated value of xT deviates from 
the experimental values. 
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V. CONCLUSION 



The calculated spin lattice relaxation rate of a local probe, with gyromagnetic ratio 7, 
in the isotropic HSM follows Eq. (|13|). This result remains valid assuming an isotropic 
interaction between the probe's spin and the molecular spin, and assuming a finite lifetime 
T~s,m for the spin state | S, m > at all temperatures and fields. Eq. (|1^) indicates that the 
spin lattice relaxation rate in these molecules is independent of the probe's gyromagnetic 
ratio at high magnetic fields, and therefore the measured values of spin lattice relaxation 
by /xSR and proton-NMR scale at the same external field (at high fields), and not at equal 
Larmor frequencies. 

The calculated lifetime of the levels r int is found to vary between 7—11 nsec, which can 
be translated to a broadening of 3 — 5 mK. Therefore we expect that the interaction which 
introduces the spin dynamics in these molecules produces level broadening of the same order 
of magnitude. 

The temperature and field independent levels broadening 1 / Ti nt calculated above can be 
attributed to an interaction that does not commute with S z and induces transitions between 
the different m states. This interaction can be dipolar between neighboring molecules, hy- 
perfine between molecular and nuclear spins, crystal field higher order terms, etc. However, 
the striking fact is that the lifetime is similar in all three molecules, indicating that it does 
not depend strongly on the spin of the molecule or the coupling J between ions inside the 
molecule, which varies greatly between the three molecules. 

This indicates that the weak dependence of the broadening on the spin value cannot be 
explained by interactions which are quadratic in S or have higher 5* dependence. This rules 
out dipolar interactions between neighboring molecules since in the three compounds the 
nearest neighbor distance is ~ 15 A. Similarly, crystal field terms which are allowed by the 
octahedral symmetry (S 2 or higheiil) are unlikely. 

The only mechanism suggested to date for level broadening of HSM, which depends weakly 
on S is the hyperfine interaction between nuclear and electronic spins. This mechanism can 
account for the finite spin lattice relaxation rate at very low temperatures!!. However, 
the values of broadening calculated above might be inaccurate due to the simplifications 
made in the calculations, but give the right order of magnitude expected from hyperfine 
interactions^. 



11 



Hyperfine interactions in anisotropic high spin molecules were studied recentlyEj'Ej, and 
their effect on QTM is becoming clearer!!!. We believe that this interaction also governs 
the spin dynamics of the isotropic molecules at very low temperatures (T < 3 K), while at 
high temperature (T > 10 K) the molecular spin dynamics are governed by spin-phonon 
interactions. 
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